Mapping the spatial distribution of entanglement in optical lattices 
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In this work we study the entangled states that can be created in bipartite two-dimensional optical 
lattices loaded with ultracold atoms. We show that using only two sets of measurements it is possible 
to compute a set of entanglement witness operators distributed over arbitrarily large regions of the 
lattice, and use these witnesses to produce two-dimensional plots of the entanglement content of 
these states. We also discuss the influence of noise on the states and on the witnesses, as well as 
the relation to ongoing experiments. 
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I. INTRODUCTION 

The quantum engineering of useful many-body states 
and the characterization of their entanglement proper- 
ties are two of the most challenging topics in Quantum 
Information Science, both theoretically and experimen- 
tally. In the laboratory, the creation of entangled states 
has been addressed in two ways. The first one starts from 
the control of individual quantum systems, let it be pho- 
tons [THl] , neutral atoms [5] [6] or ions [71 [8] , and aims 
at the creation of larger and larger many-body states. 
The second one consists of taking large numbers of these 
components (e.g., 10 3 — 10 6 atoms) and studying collec- 
tive degrees of freedom [9]. It thus seems that one has 
to make a compromise between having large entangled 
states or having a fine-grane knowledge of the properties 
of the state. 

In this work we show that there is an intermediate 
approach, by which it is possible to gain local informa- 
tion about a very large entangled state. More precisely, 
we introduce a family of operators that allow for obtain- 
ing lower bounds on the fidelity or detecting multipar- 
tite entanglement in regions of a two-dimensional graph 
state. The entanglement witnesses [T0HT4] are optimized 
for setups with ultracold atoms in 2D bipartite lattices, 
in which one now has access to the state of individual 
atoms [TSJ [16] . Remarkably, our witnesses only require 
the simultaneous measurement of all atoms, but with a 
postprocessing of the measurement statistics it provides 
a map of the quality and multipartite entanglement of 
the many-body state. 

This paper is structured as follows: in Sec. [TT] we 
review the experimental techniques available to create 
graph states in optical lattices. We present observables 
that characterize the state and act as entanglement wit- 
nesses in section III, in which we also face the main dif- 
ficulties associated to this method. We find that even 




FIG. 1: (Color online) (a) We work with two species of atoms, 
Cs and Li, trapped in two independent triangular sublattices, 
which together form a honeycomb lattice, (b) A graph state 
can be generated moving one of the sublattices along three 
different directions, (c) In this work we analyze the proper- 
ties of localizable multipartite entanglement in small sets of 
2 (a), 4 (/3), 6 (7), or more spins. Each of these regions, Q, 
is connected with its boundary, $Q, by two-qubit unitaries. 



under decoherence states useful for quantum computa- 
tion can be found, and analyse simple observables that 



bound the fidelity of the state. Finally, in Sec. IV, we 



perform numerical simulations of a cluster state subject 
to different noise sources, demonstrating that the entan- 
glement witness is capable of detecting those errors. 



II. EXPERIMENTAL GENERATION OF 
STABILIZER STATES 



The original method for creating graph states with 
neutral atoms [17] was based on filling state-dependent 
lattices with one atomic species and using controlled col- 
lisions p~8j [19]. In contrast, we will develop our ideas 
building on the experimental setup from Ref. [16] [20] , 
which traps two different species of atoms in two coex- 
isting optical lattices, one of which can be moved. This 
setup combines a diffraction mask with a powerful mi- 
croscope objective, which projects two similar triangu- 
lar lattice patterns on its focal plane. Using two light 
beams with different frequencies, the experiment may 
trap lithum and cesium atoms in two independent lat- 
tices that can be moved at will along the plane that 
confines the atoms. As shown in Fig. [ljt, we can con- 
template the Cs and Li arrangements as the triangular 
sublattices of a larger honeycomb lattice where each Cs 
atom is surrounded by three Li atoms, and vice versa; 
and each atom acts as one qubit. Since our lattice is 
bipartite by construction entanglement can be created 
using a small number of steps, equal to the coordination 
number of the full lattice. Continuing with this exam- 
ple, one has to move one sublattice three times so that 
each Cs atom approaches each of its neighboring Lithium 
atoms [Fig. [lb], suffering a controlled collision [18] or 
an engineered interaction [16 . A fundamental difference 
with previous setups [T7UT9] is that the sublattice now 
moves as a whole, regardless of the internal states of the 
atoms. If the lattices are very deep and the atom- atom 
interaction is strong enough, this can be done with great 
precision. 

To fix ideas we will assume that the entangling opera- 
tion between atoms in different sublattices is a CZ gate, 
Ucz = exp (— ^f crcs^Li) • After three parallel sets of op- 
erations, beginning with a product state, (|0) + ll))®^, 
we will arrive at a graph state 



iG )~n n 4t(io>+ii» 



^®N A +N B 



(i) 



ieAjeAf(i) 



where A and B denote the Cs and Li sublattices and 
M{i) C B is the set of nearest neighbors to the poten- 
tial well i. Note that if instead of using the control-phase 
one implements a control-NOT, £7cnot = (1 + cr£i) — 
(1 — cr^) <jg s , where the Cs absorbs the parity of its 
neighbors, we obtain what we call a "parity" multipartite 
entangled state 
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(2) 



III. ENTANGLEMENT WITNESSES 

A. Global fidelity of graph states 

All the states that we can create using the previous 
operations belong to the family of stabilizer states. In 
both cases we have a complete set of Na + Nb local 
observables, the stabilizing operators, gi, that may take 
values { — 1,+1}, and for which the states Go and Pq 
are eigenstates with eigenvalue +1 on all sites [28 . For 
instance, in the case of the graph state we have 



9i\G ) = +l\G ), WieAuB. 



(3) 



with the stabilizing operators gi = af YljeAfU) a j- ^- n & en ~ 
eral, given a set of lattice sites, Q, we can construct a 
projector onto a stabilizer state containing those sites 



^ = II*( 



l+9i)- 



(4) 



ien ' 



In theory we can use this projector to compute the fidelity 
of our experimentally realized state, p, which is probably 
mixed, with respect to the objective Gq or Pq, 



Faub =tr(P AuB p), 



(5) 



ieAjeAf(i) 



where the region under study now encloses the A and B 
sublattices. However in practice this is already impossi- 
ble for a few qubits, since the evaluation of Fq requires 
us to measure 2 Na+Nb different observables coming from 
all possible products of the gi operators. The difficulty 
of this task seems to be tantamount to performing a full 
tomography of the mixed state p. 



B. Localizable fidelity 

Instead of following this very complicated route, we 
will focus on two simpler questions, which are intimately 
related: (i) a notion of local fidelity to the stabilizer state 
and (ii) the existence and detection of genuine multipar- 
tite entanglement [21j[22] in the lattice. In both cases we 
can extract a number, let it be a fidelity or the expec- 
tation value of an entanglement witness, F(i) or W(i), 
which is distributed over the 2D lattice of sites. With 
those numbers we can study the distribution of entangle- 
ment and how much our state has been affected by noise 
or decoherence. 

Our notion of "localizable fidelity" builds on the fact 
that given a simply connected set of sites, £7, and a per- 
fect graph state, \Gq) , we can extract another perfect 
graph state in that region. One way to achieve this is 
measuring the boundary qubits, dfl, [See Fig. [lp] and, 
depending on the outcome of those measurements, per- 
forming phase gates on the qubits that were immediately 
connected to them. An alternative but completely equiv- 
alent way is to disentangle the boundary with the same 



two-qubit unitaries we used to build the state 



IV. SIMULATIONS 



pn = tr 



n n u^fpAuB 

iieonjeJV(i) 



(6) 



The most important idea is that this procedure still can 
be applied if the initial state of the atomic ensemble is 
mixed, Paub, due to decoherence. In this case the fidelity 
of the final state is related to the same observable that 
we found before, that is 



F Q = (G Q \ pa \G n ) = tr (PqPaub) , 



(7) 



the fidelity of the final state only depends on how close 
Paub is to the eigenstates of the stabilizing operators 
that cover the region and the boundary, QUdQ. The final 
observation is that the fidelity Fa gives us not only local 
information about how close our state is to the graph 
state, but it is also a witness for genuine multipartite 
entanglement in that region, Wa = ht — Pa [23] . 



C. Optimized witnesses 

However, even if (W) < detects entanglement, the 
evaluation of this quantity seems to require a number of 
measurements that increases exponentially with the num- 
ber of qubits. We thus need another ingredient, which 
is obtained by writing the fidelity as a product of two 
operators constructed from stabilizing operators corre- 
sponding to different sublattices, Pa = PanAPnnB, and 
introducing a new operator [24] 



Pn = PnnA + PanB - 1 < Pa, 



(8) 



This equation can be readily verified in the basis that 
diagonalizes both PanA and PanB, where the eigenvalues 
of the projectors can only be or +1. 

This observable provides a lower bound for the fidelity 

Fa > (Pa), (9) 

and can be used to construct an entanglement witness 



W n 



1 



-1 



(10) 



The advantage is that now the quantities (PanA) and 
(PnnB) can be extracted from just two settings of mea- 
surements. In particular, for the graph state one such 
expectation values 



(Pqha) 




n 

jeN(i) 




(11) 



is obtained measuring a x in all Cs atoms (i G O Pi A) and 
a z in the Li (j G N(i)), while the other expectation value 
is obtained with the opposite measurement basis. Note 
also that postprocessing the same set of measurement 
results we can compute the values (Pa) for any region Q, 
which allows us to produce two-dimensional displays of 
the distribution of localizable fidelity or the multipartite 
entanglement witness. 



This section features a numerical simulation of a real- 
isation of our method in an experiment taking possible 
practical sources of error into account. We have stud- 
ied the degradation of the expectation value of the wit- 
ness Wa given in Eq. (10). In general it is not possi- 
ble to compute easily the change of (Pa), but we will 
take advantage of the facts that the witness is the sum 
of two functions of stabilizing operators corresponding 
to different sublattices and that for our sources of noise 
these expectation values have simple expressions, such 
as (PanA,B) = YlieQnA^i 1 + (#*)) A in which only the 
expectation values of isolated stabilizers appear, (gi). As 
explained in Appenndix [A] we have considered various 
types of noise [25 using the quantum channel formalism 
to compute the changes in (gi). 

(i) Dephasing, which is due to fluctuations in the en- 
ergy levels of the atoms due to external fields, 
€j(p) — J dOj exp(— i<j z 9j)pex.-p(i<j z 9j)p(9j). This 
map is repeated on all sites, with site de- 
pendent uniformly distributed random phases in 
[— 6j,6j], degrading the stabilizing operator (g$) — >• 
<ft>n ie nsin(2 ei )/2 ei . ' 

(ii) Imperfections in the gates that entangle pairs of 



T U,k) 



tU*) § 



sites, Uqz* — » Uqz exp(i9jkCTjCr%,), where 9jk are 
again random variables, uniformly distributed in 
[—tjk, e jk]- This error introduces a factor in the ex- 
pectation value of the stabilizing operators, (gi) — >• 
(9i) YljeN(i) K 1 + sin(2e^)/2e^), plus other terms 



that do not contribute to the witness (10) 



(iii) Atom loss, which introduces a new state in the lat- 
tice, the hole \h). In practice, it can be described 
bycAL(p) = (l-p)p + p|0)(0|. 

(iv) Spontaneous emission €se(p) — (1 — P)p J rp\fy (0| . 

(v) The completely depolarizing channel cdp(p) = (1 — 



P)p- 



1. 



The last three sources of error have the same effect, 
(gi) -+(l-p)(gi). 

With these types of noise and decoherence we stud- 
ied the evolution of our witness operators and the overall 
description of a potential experiment using them. The 
results are shown in Figs. [2^i-c, where we plot the val- 
ues of W a , W/3 and W 1 , interpolated using smooth func- 
tions that are centered and cover the affected regions, 
a, /3 and 7, of two, four and six qubits. The result is a 
two-dimensional map of the entanglement content, where 
the value of the witness is color-coded either (a) on the 
link between two atoms, a, (b) on the central atom and 
the star surrounding it, /?, or (c) on the center of the six- 
atom plaquette, 7. In these particular plots we have com- 
bined all sources of decoherence, making some of them 
more relevant in different regions of the lattice. We have 
introduced a region of atoms subject to strong dephasing 
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FIG. 2: (Color online) Two-dimensional distribution of the entanglement witnesses for (a) two, (b) four and (c) six particle 
arrangements, a,/3 and 7 from Fig. [I]:, respectively. The value of the witness is color coded on the (a) links, (b) atoms or (c) 
center of the plaquette. A a negative value of the witness (red) denotes existence of bi- or multi-partite entanglement. All 
pictures present the same defects, consisting on two empty sites (triangles), atoms subject to strong dephasing (circle), and an 
increase of phase gate errors towards the edges of the trap. 



induced by a focused laser, covering the area marked by 
a circle. We have also emptied two sites, surrounded by 
a triangle. These empty holes are numerically equivalent 
to having spontaneous emission with 100% probability. 
Finally, we have assumed that the phase gate is 100% 
accurate in the center of the lattice and acquires a 10% 
error at the boundary of the lattice. 

We already appreciate interesting features in these sim- 
ple simulations. The first one is that bipartite entangle- 
ment is less affected by noise than multiqubit arrange- 
ments. While we can reconstruct a Bell state close to 
the boundary with an 80% error, the four- and six-qubit 
states only have an appreciable value of the witness when 
the CZ gate is above 90% fidelity. The second feature is 
that the effect of local errors remains local. The sites, 
bonds and plaquettes that share one or more qubits with 
the regions affected by atom loss or strong decoherence 
(circle and triangle in the plot) have positive values of the 
witness and do not have significant entanglement. How- 
ever, one site or plaquette away from the region of in- 
fluence of those defects, the witnesses recover their large 
negative value. 



V. CONCLUSIONS 

In summary, we have presented a simple scheme for de- 
tecting bipartite and multipartite entanglement in two- 
dimensional lattices with ultracold atoms. The present 
study admits a straightforward generalization not only to 
other bipartite lattice setups, such as square lattices, but 
also to other interaction schemes (C/cnot), or to displac- 
ing each Cs atom not three, but one or two times. First of 
all, if the Cs atoms move along two directions, the result 
is an array of linear cluster states, with an entanglement 
witness that is a generalization of the previous ones, and 



that again relies only on two-measurement settings [23] . 
If instead we move each Cs atom only once, then the 
Cs-Li interact in pairs forming a macroscopic number of 
disconnected two-qubit singlets. In this case we do not 
need a witness but can rather compute the expectation 
value of the projector 



1 



(1- 






(12) 



using three experimental settings. 

We must remark that our scheme only uses the fact 
that the lattice is bipartite and that it is possible to si- 
multaneously measure the state of all lattice sites in each 
sublattice independently. In particular, while we have 
focused on a two-species setup [16] [20] , exactly the same 
protocols and measurement schemes can be used using 
the state-dependent optical lattices in previous experi- 
ments [18], combined with the new optics that allows 
imaging individual lattice sites [26]. The only difference 
is that since we do not have different atoms on different 
sublattices, the measurement protocol has to be preceded 
by a global and local rotation of one sublattice to change 
its measurement basis. This is not too complicated and 
can be done using two counterpropagating laser beams in 
an optical lattice configuration, such that their maxima 
of intensity concide with just one sublattice. 

Our proposal represents the first experimentally real- 
izable scheme for mapping out the entanglement distri- 
bution and fidelity of a very large many-body correlated 
state. It also opens the path for the experimental detec- 
tion of very large cluster states, a task which so far was 
not achievable using ultracold atoms in optical lattices, 
but which becomes possible for ongoing experiments us- 
ing two species of atoms and holographically generated 
trapping potentials [16]. In particular, we want to remark 
that the family of graph states in honeycomb lattices is a 



universal resource for measurement-based quantum com- 
putation and that our scheme can be used to isolate re- 
gions of high fidelity in such resources. 
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1. Dephasing 

In this noise source we have an average over random 
phases 

e(cr x )= J V^Ve^p(0)#. 

If the distribution p((j)) is symmetric, then 

e(a x ) = [ [cos(20) + i sm(2(/>)a z ] a x p((j))d(j) = (1 - e { )a x , 



with some error factor e^. Since the cr z operators are not 
affected, it is legitimate to say that the map induces a 
global change in the expected value gi — >■ (1 — e^)<^. 



Appendix A: Positive maps and Noise sources 

Any physical operation on a quantum state must be a 
trace preserving positive map, which maps density ma- 
trices into density matrices. Furthermore, such operators 
admit a unique decomposition using a set of operators 



<p) = J2 Ak ? A l 



with the property e(t) = 1. 

This description admits a generalization to expecta- 
tion values. In other words, we also have a positive 
map description in the Heisenberg picture, where opera- 
tors/observables and not states are changed. Using the 
definition (6) = tr{Op}, the change in the expectation 
value can be expressed as 

(e)e( P )=tr{GY,A kP Ai} = (e(0)) 



where e(e) = Zk A { e A k- 

We now want to estimate the effect of different positive 
maps on our entanglement witnesses. We will first focus 
on local error sources. It is important to observe that the 
expected values we want to calculate are of the general 
form 

P = f(<7?GA> 4b)' a ' P e i X ' *}> 

that is, they are functions of the same observables on each 
sublattice. This means that under local error sources the 
following relation applies: 



P' = f(iKeA),i(^e B ))- 



Therefore it suffices to compute how the operators change 
under the different local error sources. In each case, the 
decoherence channel will change the effective value of the 
stabilizer expectation value. 



2. Particle loss 

This positive map has the form 

e(p) = (l-p)p + p|0)(0|, 
which we can also write in Kraus form using the operators 

A = (1 -p)l,Ai =p|0)<0|,A 2 =p|0)(l|. 
This means that the operators transform as 
e(S) = (l-p)6+p(0|6|0)l. 
Thus the stabilizer operators are modified 

g^(l-p) N a x [] *] + <£, 

jeN(i) 

where N is the number of qubits in the stabilizer operator 
(4 in our case for the honeycomb lattice) and the gf- 
contain terms that are going to vanish because they can 
be written in the form g t af HjeN>(i) a j or 9% IljeN>(i) a j 
with N'(i) C iV(i), so that their expectation values are 
zero. 

3. Errors in the gates 

We can proceed similarly though some subtleties are to 
be taken into account. First of all we realize that instead 
of transforming the state we can transform the stabilizer 
operators which appear in the expectation value 



9j ~> e 



i 52keN(j) e J kC7 j (7 k q .g+* ^ZkeN(j) e jk (J j (7 k ^ 



It can be seen that this is equivalent to performing the 
same transformation only on the a x operator 

(J* -+ Yl [cos(2e jk ) +i8m(2e jk )a?<T%] a x , 



Note that since we only have cr x operators in one sublat- 
tice and cr z on the other, the phases that we have here 
are uncorrected among different <j x operators. Further- 
more, any term that contains a <j z operator vanishes once 



we take expectation values, which means that we can re- 
place crj -)► (1 - e^crj, where €j = H k J e jk p jk (e jk )de jk . 
This shows that the outcome is a global reduction of the 
stabilizer expectation value. 



[i 

[2: 

[3; 
[4: 
[5: 

[7: 
[s; 

[9 

[10: 

[11 
[12; 

[13: 

[14 



J. Pan, D. Bouwmeester, M. Daniell, H. Weinfurter, and 
A. Zeilinger, Nature 403, 515 (2000). 
M. Bourennane, M. Eibl, C. Kurtsiefer, S. Gaertner, 
H. Weinfurter, O. Giihne, P. Hyllus, D. BruB, M. Lewen- 
stein, and A. Sanpera, Phys. Rev. Lett. 92, 087902 
(2004). 

P. Walther, K. J. Resch, T. Rudolph, E. Schenck, H. We- 
infurter, V. Vedral, M. Aspelmeyer, and A. Zeilinger, Na- 
ture 434, 169 (2005). 

C. Lu, X. Zhou, O. Giihne, W. Gao, J. Zhang, Z. Yuan, 
A. Goebel, T. Yang, and J. Pan, Nature Phys. 3, 91 
(2007). 

A. Widera, O. Mandel, M. Greiner, S. Kreim, T. W. 
Hansen, and I. Bloch, Phys. Rev. Lett. 92, 160406 (2004). 
M. Anderlini, P. J. Lee, B. L. Brown, J. Sebby-Strabley, 
W. D. Phillips, and J. V. Porto, Nature 448, 452 (2007). 
H. Haffner, W. Hansel, C. F. Roos, J. Benhelm, D. Chek- 
Al-Kar, M. Chwalla, T. Korber, U. D. Rapol, M. Riebe, 
P. O. Schmidt, et al., Nature 438, 643 (2005). 

D. Leibfried, E. Knill, S. Seidelin, J. Britton, R. B. 
Blakestad, J. Chiaverini, D. B. Hume, W. M. Itano, J. D. 
Jost, C. Langer, et al., Nature 438, 639 (2005). 

B. Julsgaard, A. Kozhekin, and E. S. Polzik, Nature 413, 
400 (2001). 

M. Horodecki, P. Horodecki, and R. Horodecki, Phys. 
Lett. A 223, 1 (1996). 
B. Terhal, Phys. Lett. A 271, 
M. Lewenstein, B. Kraus, J. 
Phys. Rev. A 62, 52310 (2000). 

J. Cirac, P. Horodecki, F. Hulpke, B. Kraus, M. Lewen- 
stein, A. Sanpera, and D. BruB, J. Mod. Opt. 49, 1399 
(2002). 
O. Giihne, P. Hyllus, D. BruB, A. Ekert, M. Lewenstein, 



319 (2000). 

Cirac, and P. Horodecki, 



[15 
[16 
[17 

[is; 

[19 
[20; 

[21 

[22 

[23; 

[24 
[25' 

[26 

[27 
[28 



C. Macchiavello, and A. Sanpera, Phys. Rev. A 66, 62305 
(2002), ISSN 1094-1622. 

W. S. Bakr, J. I. Gillen, A. Peng, S. Foiling, and 

M. Greiner, Nature 462, 74 (2009). 

K.-A. B. Soderberg, N. Gemelke, and C. Chin, New J. 

Phys. 11, 055022 (2009). 

H. J. Briegel and R. Raussendorf, Phys. Rev. Lett. 86, 

910 (2001). 

O. Mandel, M. Greiner, A. Widera, T. Rom, T. W. 

Hansen, and I. Bloch, Nature 425, 937 (2003). 

D. Jaksch, H.-J. Briegel, J. I. Cirac, C. W. Gardiner, and 
P. Zoller, Phys. Rev. Lett. 82, 1975 (1999). 

A. Klinger, S. Degenkolb, N. Gemelke, K. Brickman 
Soderberg, and C. Chin, Rev. Sci. Instrum. 81, 013109 
(2010), 0909.2475. 

C. A. Sackett, D. Kielpinski, B. E. King, C. Langer, 
V. Meyer, C. J. Myatt, M. Rowe, Q. A. Turchette, W. M. 
Itano, D. J. Wineland, et al., Nature (London) 404, 256 
(2000). 

A. Acin, D. BruB, M. Lewenstein, and A. Sanpera, Phys. 
Rev. Lett. 87, 040401 (2001), arXiv:quant-ph/0103025. 
G. Toth and O. Giihne, Phys. Rev. Lett. 94, 060501 
(2005). 

G. Toth and O. Giihne, Phys. Rev. A 72, 022340 (2005). 
G. Benenti, G. Csati, and G. Strini, Principles of Quan- 
tum Computation and Information, vol. II (World Scien- 
tific, Singapore, 2007). 

W. S. Bakr, J. I. Gillen, A. Peng, S. Foiling, and 
M. Greiner, Nature 462, 74 (2009). 

D. Gottesman, Phys. Rev. A 57, 127 (1998). 

The gi operators are the generators of the so called sta- 
bilizer group |27| . 



